QED of lossy cavities: operator and quantum- state input-output relations 
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Within the framework of exact quantization of the electromagnetic field in dispersing and ab- 
sorbing media the input-output problem of a high-Q cavity is studied, with special emphasis on the 
absorption losses in the coupling mirror. As expected, the cavity modes are found to obey quantum 
Langevin equations, which could be also obtained from quantum noise theories, by appropriately 
coupling the cavity modes to dissipative systems, including the effect of the mirror-assisted absorp- 
tion losses. On the contrary, the operator input-output relations obtained in this way would be 
incomplete in general, as the exact calculation shows. On the basis of the operator input-output 
relations the problem of extracting the quantum state of an initially excited cavity mode is studied 
and input-output relations for the s-parameterized phase-space function are derived, with special 
emphasis on the relation between the Wigner functions of the quantum states of the outgoing field 
and the cavity field. 

PACS numbers: 42.50.Dv, 42.50.Lc, 05.30.-d 



I. INTRODUCTION 

The use of atoms interacting with light has been very 
promising in handling information storage, communica- 
tion, and computation 0, 0> EJ 01 • I n fact, optical sys- 
tems do not only allow the observation of fundamental 
quantum effects j(| 0, H, , but they can be also used to 
implement quantum networks with photons, which may 
be regarded as representing the best qubit carriers for 
fast and long-distance quantum communication |lOUll| . 
In optical systems resonatorlike devices — referred to as 
cavities in the following — are indispensable elements. In 
particular, high-Q cavities have been well known to offer 
a number of possibilities to engineer nonclassical states 
of light [UGll- 

Since high-Q cavities feature well-pronounced line 
spectra of the electromagnetic field, which renders it pos- 
sible to control the atom-field interaction to a high de- 
gree, they are best suited for the generation of quantum 
states on demand. For example, in Ref. 01 the creation 
of arbitrary superposition of Fock states of light inside a 
cavity by means of controlling the time sequence of the 
amplitudes and phases of the atom-field interactions of 
the quantized cavity field and the external driving fields 
with a trapped atom is considered. Another proposal 
to generate superposition Fock states or coherent states 
in a cavity is to exploit adiabatic interaction of the cav- 
ity field with an atomic system by achieving the transfer 
of ground-state Zeeman coherence onto the cavity-mode 
field 0]. It is worth noting that the idea of the gener- 
ation of a bit-stream of single photons on demand in an 
optical cavity is based on the concept of adiabatic pas- 
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sage 0, 0, 0. In the specific context of microwave 
cavities, a scheme proposed in for the generation of 
photon number states on demand via it pulse interaction 
of single two-level atoms passing through a cavity, has 
been experimentally realized d|20t l2l|. 

To measure the quantum state of a cavity field, vari- 
ous schemes have been considered. As has been demon- 
strated experimentally [2^, |23| , the quantum state of a 
microwave cavity field can be reconstructed by employ- 
ing the dispersive interaction of a single circular Rydberg 
atom with the cavity field [24| ■ To reconstruct the quan- 
tum state of a cavity field by measuring the field escap- 
ing from the cavity, the proposal has been made to use 
pulsed homodyne detection and an operational definition 
of the Wigner function in terms of appropriately chosen 
collective mode operators plf . 

However, to further use intracavity-generated photonic 
quantum states in applications such as quantum networks 
connecting distant quantum processors and memories, 
the issue of extraction from the cavities of the quantum 
states needs to be considered very carefully with resp ect 
to quantum decoherence. In the scheme in Refs. |2(| l27| 
qubits that are stored in the internal states of cold atoms, 
located in the antinodes of a standing wave of a high-Q 
optical cavity are mapped onto photon number states, 
which play the role of a communication channel by leak- 
ing out of the cavity and being caught in a second cav- 
ity. Further, schemes for the generation of entangled 
states of individual atoms held in distant cavities have 
been considered (see, e.g., Refs. (2£| H|| |3(j ) ■ In all the 
schemes the absorption losses unavoidably occurring in 
the processes of exit and entrance of the photons through 
the coupling mirrors are typically disregarded. However, 
even if from the point of view of classical optics these 
unwanted losses are very small, so that they effectively 
do not influence classical light, they can lead to a drastic 
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degradation of nonclassical light features indispensable 
for quantum communication (see, e.g., Ref. HJl). 

Roughly speaking, there have been two routes of treat- 
ing the input-output problem of a leaky cavity. In the 
first — the quantum stochastic approach to the problem — 
standard Markovian damping theory is employed, where 
the dynamical system is identified with a chosen mode 
of the perfect cavity ( Q— >oo), the dissipative system is 
identified with the continuum of modes outside the cav- 
ity, and a bilinear coupling e nerg y between the modes of 
the two systems is assumed 32]. In this way, the cav- 
ity mode is found to obey a quantum Langevin equation, 
and operator input-output relations can be derived. The 
theory can be used, e.g., to relate correlation functions 
of the outgoing field to correlation functions of the cavity 
field and the incoming field or to describe the cou- 
pling of modes of two cavities through their respective 
input and output ports |34l l35| . 

In the second route — the quantum field theoretical ap- 
proach to the problem — the calculations are based on 
Maxwell's equations and exact quantization of the elec- 
tromagnetic field in the presence of nonabsorbing cav- 
ity walls described in terms of appropriately chosen real 
permittivities [IE S3 • Having established the equiv- 
alence of the two routes, thereby constructing the in- 
teraction energy between the cavity field and the outer 
field, one may try to include unwanted absorption losses 
in the theory by allowing for further dissipative systems 
in such a way that appropriately chosen interaction ener- 
gies between them and the cavity modes are added to the 
Hamiltonian used in the quantum stochastic approach 
[38l l39l | . As we will show within the frame of exact quan- 
tum electrodynamics in dispersing and absorbing media, 
this simple concept, though leading to the correct quan- 
tum Langevin equations for the cavity field, does not lead 
to the correct operator input-output relations in general, 
because the absorption losses in the coupling mirrors are 
not properly taken into account. 

If the operator input-output relations are known, the 
correlation functions of the outgoing field can be ex- 
pressed in terms of correlation functions of the cavity 
field and the incoming field. In this context the question 
of the calculation of outgoing-field quantum state as a 
whole arises. Starting with the correct operator input- 
output relations that include both wanted and unwanted 
losses, we will calculate the quantum state of the pulse- 
like field escaping from a high-Q cavity, assuming that 
the quantum state of the corresponding cavity field at 
some initial time is known. 

The paper is organized as follows. In Sec. ITT1 some basic 
equations are given and the cavity model is introduced. 
The intracavity field and the outgoing field, including the 
operator input-output relations, are studied in Sees. IIIII 
and IIVI respectively. The problem of extraction of an 
initially prepared cavity-quantum state is considered in 
Sec.0 and a summary and concluding remarks are given 
in Sec. IVII Some derivations are given in appendices. 



II. PRELIMINARIES 



Quantization scheme 



Let us consider N atoms (with the Ath atom being 
at position ta) that in electric dipole approximation in- 
teract with the electromagnetic field in the presence of 
linear dielectric media of spatially varying and frequency- 
dependent complex permittivity 



e(r, bj) = s'(r, lu) + is"(r, u>). 



(1) 



Note that due to causality the real and imaginary parts 
e'(r, u) and e"(r, lj), respectively, are uniquely related to 
each other through the Kramers-Kronig relations. Fol- 
lowing the approach to quantization of the macroscopic 
Maxwell field as given in Refs. [U EJ IS |H , we may 
write the multipolar-coupling Hamiltonian in the form of 

M 



H = H, 



field 



+ H: 
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Here. 



fffieid = J d3r J o dwfiwft(r,w)-f(r,( 



(2) 



(3) 



is the Hamiltonian of the system composed of the electro- 
magnetic field and the medium, including a reservoir nec- 
essarily associated with material absorption, with f (r, uj) 
[and F(r,w)] being bosonic fields that play the role of 
the dynamical variables of the composed system, 

[/ M (r,uO,£,(rV)] = 5^5{lo-J)5^{v-v>), (4) 



[/ M (r,w),/ M ,(r>')] =0 



(5) 



(the Greek letters label the Cartesian components). Fur- 
ther, 



-Hatom — ^ ^ f^AkSAkk 
A k 

is the atomic Hamiltonian and 



Hint = -^2d A -E(r A ) 



(6) 



(7) 



is the (multipolar-)interaction energy, where 

S Ak 'k = \k') AA {k\ (8) 
are the flip operators of the Ath atom, 

dA = dAkk'SAkk' (9) 



kk' 



is its electric dipole moment (dAkk> = A{k\dA\k') a), and 
E(r) is the medium-assisted electric field, which ex- 
pressed in terms of f(r,cj) [and f^(r,o;)] reads 

/>oo 

E(r) = / dwE(r,w)+H.c, (10) 
Jo 
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E(r,w) = i^ d — lo 



0+ 



x / d 3 r'y/£"(r', u) G(r, r', w) • f (r', w) (11) 



where the classical Green tensor G(r, r',oj), which also 
corresponds to the quantum field-theoretical retarded 
Green tensor (see, e.g., Ref. 46] ), is the solution to the 
equation 

(.2 

V x V x G(r, r', - — e(r, w)G(r, r', w) = <5 (3) (r - r') 
together with the boundary condition 



G(r, r', ui) — > for |r — r'| — > oo 



(13) 



It is not difficult to see, that in the Heiscnberg picture 
f(r, u>,t) obeys the equation of motion 

f(r,w,t) = —\f(r,u,t),H] = -iwf(r,u,t) 
in L J 

+ /WW^| y^^) ^d A (t) ■ G*(r A ,r,u), (14) 

A 

which after formal integration leads to 

f(r,w,t) = ffree(r,w,t) + f s (r,w,t), (15) 

where 

f frcc (r, u, t) = e-^(*- t ')f free (r, w, (16) 

and 

f s (r,w,i) = hquj 2 \ /l^ y/e"(r,u) 

V /17T 

x^ /'dt'e(t-OdA(t')-G*(rA,r, W ) e -^ t -*''. (17) 

Substitution of Eq. (JT^J together with Eqs. and fP7ll 
into Eq. (|ll|l yields the corresponding source-quantity 
representation of E(r, ui, t). 

B. Cavity model 

For the sake of transparency, let us consider a one- 
dimensional cavity modeled by a planar dielectric 4-layer 
system (Fig. EJ. In particular, the layers j = and j = 2, 
respectively, are assumed to correspond to perfectly and 
fractionally reflecting mirrors which confine the cavity 
(layer j = 1). In what follows we use, with respect to z, 
shifted coordinate systems such that < z < I for j = l, 
< z < d for j — 2, and < z < oo for j — 3. Applying 
the one-dimensional version of Eq. (|llf> together with 
Eqs. 115|) - (|17fl to the field in the jth layer of permittivity 
£j(ui) (j = 1, 2, 3), we may write 



m (z, w, t) = ti£jz, u>, t) + E^'(z, u>, t), (18) 



U), 
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FIG. 1: Scheme of the cavity. The fractionally transparent 
mirror [region (2)] is modeled by a dielectric plate. The active 
sources are in region 1, which can also contain some medium. 



with 



E flcc (z,uj,t) = iujuq 



£ I dz'Gto')(z,z', W )j;2(z',c,t) (19) 
_ r =i J W] 



and 



E / d*'e(*-t') 



x e-^-^dAit^lmG^izAiZ^) (20) 

(.4, mirror area), where the integral relation ()D2fl has 
been employed. Here, [j 1 ] indicates integration over the 
j'th layer, the abbreviating notation 



^{z,v,t)=w^s'j{w)fg e (z.^l) (21) 



is used, and it is assumed that the active atomic sources 
are localized inside the cavity. The (nonlocal part of the) 
Green function reads |47| 

G^ j '\z,z',u) 

= li[S^ > (z,uj)E"'S^ < (z\uj)Q(j - f) 
+ ui) w)Q(j' - j)] , (22) 

where the functions 

£ U)> (z, w) = eW**-*** + r j/3 e- ifi **- d * ) (23) 

and 

fW<(z,w) = e-Vt'+Tj/oeW'*, (24) 

respectively, represent waves of unit strength traveling 
rightward and leftward in the jth layer and being re- 
flected at the boundary [note that Q(j — j') means 
0(z — z') for j = j']. Further, S" is defined by 



~33 



1 to/^^i t 3/ i,e i0 i' d 



3/j' 



(25) 
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where 



and 



Da 



r j/0 r j/3 



,2il3jdj 



(26) 



f3 s = f3 s {u 3 ) = J7&)~ 



= [n>)+m?(u;)]- = $+^ (#,$'>0) (27) 

(di=l, d2=d, d 3 — 0). The quantities ij/j' = (/3j / (3j> )tj> u 
and Tj/j/ denote, respectively, the transmission and re- 
flection coefficients between the layers j' and j, which 
can be recursively determined (for recursion formulas, 
see Appendix E| . 



m»,t) = -^f m M^£Lv,»,t), (3D 



Hoc 
2m 



[3] 



dz'e^'ffjz',^), (32) 



and 



D' 2 (lo) = 1 - r 2 ir 2 3e 



2ip 2 d 



(33) 



Inspection of Eq. (|28|l shows that the function Di(ui) 
defined by Eq. I|26|) for j — 1 characterizes the spectral 
response of the cavity. In particular, its zeros determine 
the (complex) resonance frequencies Qk , 



D 1 (n k ) = l + r 13 (n k )e 2l ^ n ^ 1 =0. 



(34) 



III. CAVITY FIELD 

To further evaluate the equations given above, we first 
consider the field inside the cavity (j = 1). In order 
to make contact with the familiar standing-wave expan- 
sion in the idealized case of a lossless cavity, it is useful 
to rewrite the equations with the aim to obtain a non- 
monochromatic mode expansion that takes into account 
the finite line widths due to the wanted input-output cou- 
pling and the unwanted absorption losses that unavoid- 
ably exist in practice. 



A. Nonmonochromatic mode expansion 



We begin with the free field. From Eqs. (jT§|) and (|2*2|l it 

follows that E^ T l e (z, uj, t) can be represented in the form 
of 



E_H e {z,u),t) 
1 

D~i 



C ( H(z,u J ,t)-C ( i ) _(z,u J ,t) 



t2ie 



D' 



(7i 2) (c,i)+r 23 e 2 ^cf( W ,t) 



(28) 



where 



-M / & z >Q{ z - z >)e^'jH{z',u,t), (29) 
2ni J m - tree 



-M f Az'Q{z'-z)e^'j ( l ljz',cj,t), (30) 
2m J m - tree 



Note that when the coupling mirror is not a single plate 
but — as in practice — a multilayer system, then ri 3 (io) 
is the reflection coefficient of the multilayer system and 
Eq. (|34|l applies as well. Decomposing il k into real and 
imaginary parts according to 



(35) 



we can write the formal solution to Eq. (|34|l in the form 
of 



1 



and 



I 2 ni 



x <n 



27rfc+7r-tan" 1 4^ 



13 



-n'/lnlnal (36) 



r fc = 



c 1 



I M 2 



x <n 



2irk+ir— tan 



13 



13 



+n[ln\n 3 \\ (37) 



[ni=ni(Ofe), ria =ri 3 (Ofe)], from which w fc and T fe may 
be calculated by iteration, by starting, e.g., with the res- 
onance frequencies of the lossless cavity. 

Let s(t) be a function of time whose Fourier transform 
is given by 



(38) 



and assume that S(u>) is analytic in the lower half-plane. 
Employing the residue theorem, we may write 



W 1 2ir Di(w) 



Y — 



Q(t) ( 



(39) 

Applying Eqs. (|3"5j) and l|3"9l to the c- number func- 
tions sm(j3 1 z)D- 1 , sin(/3 1 z)r 13 e 2l ' 3li £>- 1 , and [e tf3lZ + 
ri 3 e~ ll3l ( z ~ 2[S> }D~ 1 in Eq. (|2*5|l and disregarding (irrele- 
vant high-frequency) contributions that may arise from 
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poles other than those of D 1 we may rewrite To calculate the electric free field 



E bcc (z,uj,t) as 



(1) 



(40) 



with 



E kfTCC (Z,LU,t) = I 



sin[#.(fi 



T(cj)b- m (u, t') + Mu)cx{uj, t')] (41) 



(A = cav, +, — ), where the operators c\(uj, t) are defined 
according to 



Ccav(w,t) = -a c 



/ nA /ioc 
fiochuj 1n\ 



x / dz sin(/3!z)j 
7[i] 



(42) 



c±(w,t) = a±i 



' nA 
HqcHlu 



£ / + \ 2 KI / A(3)/ ,n 

b hl (uj,t) = —j=\ — Cl'(u,t), 



Hochuj 



(43) 
(44) 



with 



"cav = «cavM 

= 2V2|ni| K sinh(2#'Z) - n" sin(2/3^)]" 
a± = a±(u>) 



(45) 



"2 6 



n' 2 smh(^'d) ± n' 2 ' sm(/3' 2 d)] 2 , (46) 



A cav (w) = -4i 



(r 23 e^ d ± 1) 



(47) 



(48) 



(49) 



It is straightforward to prove that the operators &\(u),t) 
satisfy Bose commutation relations: 

[c x (u>,t),c{,(u/Jj\ =5 xx ,6{u-u')e- i '» ( t- t '\ (50) 
[& in (a,,t),^>'.0] = .5( w - w ')e-Mt-0 (51) 

with all other commutators being zero. 



4ee(^ *) = / d ^ IfeLfo W : *) + H - C - ( 52 ) 
JO 

[cf. Eq. IjlOp]. we subdivide the ui axis into intervals 
(Afc)= [i(wfc_i +w fe ), i(w fc +cj fc+ i)] and write 



4 1 e ) e(^*) = E^ 1 fr ) ee(^*)+H.C. 



(53) 



where 



£2L(M)=/ duE^^t) (54) 

(recall that the index k is used to numerate the reso- 
nances of the cavity). Substitution of Eq. (|40|l together 
with Eq. (JUJl into Eq. (53} yields 

^fcfree ( Z > i ) = / dw ^fefree 

E / dw£^j rco (z,w,i). (55) 



For sufficiently high-Q cavities, i.e., -C Aw^, with 
Awt = — being the width of the fcth inter- 

val, the second term in Eq. I|55|l can be regarded as being 
small compared with the first one and may be omitted in 
general, leading to 

£«L(*.*)=/ Au^l cc {z,u,t). (56) 

J(A k ) 

In this approximation, Eq. I(53|) reduces to 

4oe(^ *) = E / dW ^SL^ *)• (57) 

Note that within the approximation scheme used, the 
lower (upper) limit of integration in Eq. I(57|) may be 
extended to —00 (+00). 

To determine the source field 



4 (1) M) 



dujE^\z,Lu,t) + H.c, (58) 



we start with Eq. l(2Tj|) together with Eq. (J22J. Perform- 
ing the Fourier transformation and using the resonance 
properties of the cavity response function [see Eqs. 1)34(1 - 
l)39p]. we obtain, in the same approximation that leads 
from Eq. JUJ to Eq. |57|). 



EP(z,t)=J2Eg\z,t) +H.c 



(59) 



where 



x e 



e £i(O fe )M 
«l*(*-f )^ (t /j sin^ri^ZA] + H.c. . (60) 
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Note that in Eq. ({BTijl it is assumed that e luJkt dA(t) may 
be regarded as being an effectively slowly varying quan- 
tity. Combination of E^ e (z, t) and Ei'(z,t) to the full 
intracavity field 

E^{z,t)=E<£i{z,t) + E^{z,t) (61) 
yields the nonmonochromatic mode expansion sought. 

B. Quantum Langevin equations 

To bring Eq. (gU [with E^(z,t) from Eq. JSZJ to- 
gether with Eq. l|IT)l and E^ [z, t) from Eq. together 
with Eq. PJ] in a more familiar form, we introduce the 
operators 



&kin(*) = 



1 



'27T J( Afc ) 

7=/ 

/27T J( Afc ) 



dw ca(w,^), 
dw 6in(a;,t), 



(62) 
(63) 



which, on a time scale At^> Aio k , Au> k , , obviously obey 
[recall Eqs. (|50(l and l|51|l] the commutation relations 

[c k x(t),£ly(t')] =d w S xy S(t-t% (64) 
[b ki n(t),b{, hl (t')} =5 kk ,5{t-t>). (65) 

Further, recalling Eqs. l|3T]). lt5T|) . and (JSDJl, we may 

rewrite Eq. JJHJ as (T k < Aw fc ) 



^ 1 )(^t) = ^ J E fc (^)a fe (t)+H. 



(66) 



where the standing wave mode functions are defined as 

i 

h 



E k (z) = iuj k 

and 

a k (t) -- 



e ei(uj k )lAuj k 



sm\f3i(u>k)z], (67) 



di / 0(t-f / )e~ <n * ( *~*' ) 

TUfcin(i')+5^fcA4A(0 



2ni(w fc )Z 



ft- 



(68) 



[T fc = T(w fc ), 4a^aW]- 

From Eq. (|68|l it is not difficult to see that cifc obeys 
the Langevin equation 



2m(w fc )f 



TfcS fein (t)+^Afe A CfeA(t) 

A 



(69) 



and it can be proved (see Appendix [D| that the equal- 
time commutation relation 



[6fc(*)>6fc/(*)] = <W 



(70) 



holds. 

The damping rate in the first term on the right-hand 
side of Eq. (|69|) can be decomposed as follows (see Ap- 
pendix Q : 



Tfe — 7fcrad + 7feabs , 

7ferad = 2KKW |Tfc 

7fcabs = E 7fcA 



2|m(w fc )|i 



(71) 
(72) 

(73) 



A 1 v " A 

Here, 7/crad is the radiative decay rate describing the 
transmission losses due to the input-output coupling and 
7fcabs is the (nonradiative) decay rate describing the ab- 
sorption losses inside the cavity (term proportional to 
|^4-fccav| 2 ) and inside the mirror (terms proportional to 
|Afc±| 2 ). Accordingly, the Langevin noise force as given 
by the third term on the right-hand side of Eq. H69|) 
consists of the contributions associated with the losses 
due to the input-output coupling [term proportional to 
T k b k i n (t)] and the absorption losses inside the cavity 
[term proportional to A kca-V c kca , v (t)] and inside the mirror 
[terms proportional to A k ±c k ±(t)]. 

Equation (|69|l can be regarded as a generalization of 
the results derived in Ref. J36J for a leaky cavity without 
material absorption to a realistic cavity which gives rise 
to both radiative and unwanted (nonradiative) absorp- 
tion losses. In particular, when ei(u> k ) can be regarded 
as being real, then the second term on the right-hand 
side of Eq. (|69|l is nothing but the familiar commutator 
term (ifi,) -1 [afc,-ff; n t], where 



Hint = - E E E k{ z A)dAa k + H.c. 



(74) 



Moreover, from Eq. together with Eqs. J7T}-lO it is 
seen that the effect of absorption losses on the intracavity 
field may be equivalently described within the frame of 
Markovian damping theory, with 



abs.int 



= *EE 

A 

+ H.c. 





c 


/ du 


2n\(uj)l 





A k \(u)a\c k \{uj 



(75) 

being the total interaction energy between the cavity 
modes and the dissipative systems responsible for absorp- 
tion. Thus Eq. H69J) can be also reg arded as an extension 
of the results derived in Ref. [33j within the frame of 
quantum noise theories, by adding to the Hamiltonian 
therein an interaction energy of the type (|75|) ■ Needless 
to say that also other than the dissipative channels con- 
sidered here can be included in the interaction energy. 
The unwanted losses attributed to the cavity wall that 
has been assumed to be perfectly reflecting is a typical 
example. Note that Eq. (|75|l implies that each cavity 
mode is coupled to its own dissipative systems. 
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IV. FIELD OUTSIDE THE CAVITY 

Once the cavity field is expressed in terms of non- 
monochromatic mode operators a k (t) , the question arises 
of how the outgoing field is related to it. To answer, we 
first rewrite the outgoing field using Eqs. lfTB|) - (|20)l with 
j = 3 and proceed similarly as in the case of the cavity 
field. 



A. Outgoing field 

We again begin with the free field. Inserting the Green 
tensor as given by Eq. in Eq. fli^l (j = 3) and separat- 
ing the incoming and outgoing parts propagating along 
— z and z, respectively, we may represent the outgoing 
part at z = + (cf. Fig.^l as follows (see Appendix |B|) : 



E 



(3) 



ut,frec( Z : u i *) | z =0+ 



ipil r 



D' 



+ 



D> 



where C<±(1, oj,t) is defined by Eq. lt2T?|) (for x = l), and 



2i/3 2 dn( 2 )( 

- i 3 ie 4,3l 'ci 3) (w,0 
(w, t) + r 21 ci 2) (u, tj\ +r 31 C& (oj, t), 



(76) 



C ( ±>(u>,t) and Cf } {u,t) are defined by Eqs. and 
(|32|l . respectively. Treating the term D^iuj) in the same 
way as that leading from Eq. I|28|) to Eq. I|4l)|l together 
with Eqs. (gTJ-gSJ), from Eq. we derive 



^(3) 



i?out,free( Z ' w i0| z=0 4 



X^°Ufrce(X W ^)| z=0 -M (77) 



where 



p (3) ( +\ 

£Lko\lt,hcc\ Z i W : l ) 



1 / flQCfrliJ 



2 V TrArii 2ni{fl k )l 
(w)ca(w,^') 



l fiochio 

2\ nA 



2.; 



1% 



|na| 

i/3 2 d _|_ ^ 



&fc ln (cj,£) 



C/c+(w,t) 

Cfe-(w,*) 



(78) 



As we see from this equation, there are three physically 
different contributions to the outgoing free field. The first 
(integral) term proportional to £13 represents the frac- 
tion of the cavity field transmitted through the mirror 
[cf. Eq. lHJl ]. The term proportional to r 31 b k - m (uj, t) rep- 
resents the reflected part of the incoming field, whereas 
the terms proportional to t23C k ±(io,t) describe the field 
attributed to the noise sources inside the mirror. 

Integrating Eq. I|78|) with respect to u, we obtain the 
free-field part of the outgoing electric field in the time- 
domain, 



E. 



(3) 

out, free 



z=0+ 



k 



=0+ 



-H.c. (79) 



[cf. Eqs. (|52|l and Ij53(l ]. where, within the approximation 
scheme used, 



E 



(3) 

fcout,frcc 



(-M)I 



Z = + 



.(3) 



(A fc ) 



d^^feout,free(^^*)L= 



+ 



(80) 

[cf. Eqs. |53>-(|57|>]. 

Starting from Eq. (HU (j = 3) together with Eq. J22J, 
we may rewrite the source-field part of the outgoing elec- 
tric field to obtain, in close analogy to Eqs. I|58|l and (ICiUI) . 



E. 



where 



(3) 



,(M)L 0+ =E&(z,t)\ z= 



2=0 + 



EC(m)L 0+ +h,. 



(81) 



duE ks (z,uj,t)\ z=0+ +H.c. 

(A fe ) 

u k ti3 e-xp[if3i(tt k )l] 

A 



dt'e{t-t') 



x e 



2e ei(O fe )U 
-^( i -*')^(t')sin[/3 1 (n fe )^] 



H.c. 



(82) 



Finally, combination of the free-field part and the source- 
field part yields the full outgoing field at z = Q + , 



z=0+ 



E. 



out, free ^)| z =0+ + ^s 3 '( z 7 *)| z=0 +- 

(83) 



B. Global input-output relation 

Let us restrict our attention, for simplicity, to a cavity 
in free space, i.e., n 3 — > 1, and define the operators 



£>out(w,£) = 2* 



— j-E^ ut (z^t)\^ 0+ . 



(84) 



By using the formulas given in Section flV Al it is not diffi- 
cult to see that we may rewrite the w-integrated operator 



&out(*) 



27T 



du) b ou t(u,t) 



(85) 
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&out(i) = y^fe/cout(o, 



where 



bkout(t) 



1 



'27T J (A fc ) 

with 6 fcout (w,i) being given by 



(86) 



(87) 



-'fcout 



(w,t) 



-^-^(z, W ,t)| z=0+ 



+ Ajg ( w )c fe+ (w, i) + AW ( w )c fe _ (w, i) 
+ 4 o) H^in(^,0- 



T k (w)b k in{u, f) + ^ A k \(u)c k \{u, t') 



Here, the functions A k °±(u>), R k "'(uj), and T k "'(uj) are 
defined as follows: 



?(»)/ 



t 2 3 1 ± ^2ie 



(°) 



V«4 



4 o) M = ^i. 



(89) 
(90) 
(91) 



Performing in Eq. 1(8 7() the uj integral, on extending 
again the lower (upper) limit to — oo (+oo) and recalling 
Eqs. JBSJ and JS^Jl, we see that the source term in Eq. 
and the second (integral) term in this equation sum up 
to a term proportional to the cavity- field operator a k (t) . 
Thus, from Eqs. it follows that (r fc < Au; k ) 



Jfcout 



(*) = 



2ni(uj k )l 



Tl o) a k (t)+R k o) b kin (t) 



4-4-(i) 



(92) 



If on the right-hand side in Eq. 192(1 the third term and 
the forth term, which result from the absorption losses in 
the coupling mirror, are omitted, then Eq. 19211 reduces 
to the well-known input-output relation 33, 36] for a 
leaky cavity whose losses solely result from the wanted 
radiative input-output coupling, in which case | -Hi (^fe) | 
= 1 holds. It can be shown that the operators that appear 
in Eq. (|92H obey commutation relations of the type given 
in Ref. [36j. In particular on the time scale considered, 
i.e., At^> Au> k ,Au> k , , the commutation relation 



[Stout W.SLutCO] -0 



(93) 



holds (see Appendix |EJ. 

As we know from Sec. 1111 lH the damping of the cavity 
modes due to unwanted losses can simply be described by 
introducing into the Hamiltonian an interaction energy of 
the type ((7511 and treating its effect in Markovian approx- 
imation. However, since such a simple interaction energy 
does not allow to include in the theory effects such as for 
example the influence of the coupling-mirror-assisted ab- 
sorption on the outgoing field via the incoming field, the 
last two terms in Eq. (|92l) arc missing and hence \R k \ = 1 
is set (see, e.g., Ref. [33|)- As a matter of fact, due to 
unavoidable losses it is always observed that |iZi°'| < 1 
in practice. The additional noise associated with these 
losses is just described by the last two terms in Eq. 1(92(1 . 

The input-output relation can be used, e.g., to cal- 
culate correlation functions of the outgoing field in terms 
of (in general mixed) correlation functions of the cavity 
field, the incoming field, and the dissipative channels. In 
what follows we will not consider the one or the other 
correlation function, but focus on the quantum state as 
a whole. In this connection it should be stress laid on the 
fact that Eq. ((92(1 as a global input-output relation does 
not specify the incoming and outgoing (nonmonochro- 
matic) modes as well as those of the dissipative channels 
which are really connected in the relevant frequency in- 
terval defined by the bandwidth of the cavity mode and 
hence basically carry the quantum state. 



V. QUANTUM STATE OF THE OUTGOING 
FIELD 

For the sake of transparency let us suppose that during 
the passage of atoms through the cavity the fcth cavity 
mode is prepared in some quantum state and assume 
that the preparation time is sufficiently short compared 
with the decay time 1^ , so that the two time scales are 
well distinguishable. In this case we may assume that 
at some time t$ (when the atom leaves the cavity) the 
cavity mode is prepared in a given quantum state and 
its evolution in the further course of time (i.e., for times 
t>to) can be treated as free-field evolution. To specify 
the relevant modes, it is useful not to use Eq. 1(921) but 
return to Eq. I(88|) and relate therein b kou t(uJ, t) to a k (to)- 
It can be proved (see Appendix[FJ) that on the (relevant) 
time scale At^s> Acj^T 1 Eq. 1(88(1 can be rewritten as 



2ni(u k )l 



T, 



(o) 



/27T 



t+At 



dt' e-^-^a^t') 



to 



+ A k % + (uj, i )e-^ t0 ) + A k °lc k _ (u, t )e-^-^ 



+ 4% in (^ )e- Mt -' o) . 



(94) 



Note that integration of both sides of Eq. ((94(1 with re- 
spect to to over the interval (A k ) leads to the input- 
output relation 1(92(1 (713 — ► 1). It should be mentioned 
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that, for the special case of purely radiative losses, an 
equation of the type of Eq. t9"ll could be also found from 
the quantum stochastic theory in Ref. [23, which would 
however suggest that its validity only requires the condi- 
tion At > to be satisfied. 

Substituting Eq. i jHSjl (for t > to) together with 
Eqs. and into Eq. (jMJ), we derive 



&fcout(y,i) = F k (uj,t)a k (t ) + B k (uj,t), 
where the c- number function F k (uj,t) reads 



(95) 



F k (uj,t) 



1/2 



2tt 

exp - n* k )(t + At- t )] - 1 



(96) 



and the operator B k (u,t) is a linear functional of the 
operators b kin (u,t ) and c k \(oJ,t ): 



where, for chosen t, the nonmonochromatic mode func- 
tions 4>i(uj,t) are a complete set of square integrable or- 
thonormal functions: 



(A fc ) 



du>4>i(u,t)(f>*j(LJ,t) = 5ij, 



(104) 



J2<t>i(u,M("',t) =S(u-w'). (105) 



Needless to say that the commutation relation 

KLW>ffi(*)]=*j (106) 

holds. 

Let b^ ut (t) be the operator attributed to the outgo- 
ing mode that is related to the cavity mode through the 
input-output relation (|95|l : 



F k (m,t) 



(107) 



B k (u>,t) = / da/ 
J {Ah) 



<*ki» (<*>,<*>',*) b kin (u',t Q ) 



(97) 



rj k (t) = / du\F k {u,t)\ 2 . (108) 

J{A k ) 

By using Eqs. J3SJ|, l|107|) . and (|108|l we may rewrite (for 
chosen k) Eq. (|103fl as 



Here, 



Gfci n (w, a/, t) = T^*T k Vk(uJ,Uj',t) 

G fccav (w,cj',t) = T^ o) *A£ cav w fc (a;,a/,i), 
G k± (w,w',t) =T { ° ] *Al ± v k {Lo,Lo',t) 



with 



Ice' 



-iujAt 



2tt 2nlluj-tt* k 

e iu'(t+At-t ) _ e iO*(t+At-i ) 

iw(t+At-t ) _ „«y(t+At-t ) 



(98) 
(99) 

(100) 



(101) 



V / %Wa fc (to) + ^ ) (0 ift = l, 



otherwise, 



where 



(109) 



Bf{t) = I dw&(w,t)B fc (w,t). (110) 
J(A k ) 

Inserting Eq. (|97|) in Eq. (|110|) , we may rewrite B^' (t) 



as 



^ ) (i) = v / ^^L(*)+Ev / 5f(*)4t(*)- an) 



Here the functions £^ (t) (cr = in, A) read 

C^(*)=/ do,| X «(c,t)| 2 , (112) 



where 



A. Nonmonochromatic modes 

To calculate the quantum state of the outgoing field, 
it is convenient to introduce a unitary, explicitly time- 
dependent transformation according to 



SwM) = ^*(u;,t)S» ut (f), 



(102) 



6£Lt(t) 



(A*) 



duj^i(u},t)bkovit(^,t), (103) 



Xk %,t)= du/ (113) 



The operators &jy n (i) and c k {(t) are defined by 



= / do; 



(A fc ) 



6/Hn(w,*o), (H4) 



dw 



c k \(uj,to). (115) 
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B. Phase-space functions 

Introducing in the characteristic functional 

Cfcout [fl{u)),t\ 



exp 



(A* 



dw i 9H6L ut (a;,t)-H.c. 



(116) 



the operators Vui t (t) according to Eq. 1|102|) and taking 
into account the commutation relation 1(106() we see that 
the operator exponential factorizes as 



exp 



(A k ) 



nexp[/34out(*)-H.c. 



where 



which follows from Eq. (|109|> together with the commu- 
tation relations J7UJ), ((106(1 . and 1120(1 . the further evalu- 
ation of Eq. ((122(1 is straightforward. Following Ref. jsjj 
and calculating the characteristic function C ou t(/3, t; s) = 
Cfcout (/3ij ^; s) in s order of the quantum state of the rel- 
evant outgoing field, we may express it in terms of the 
characteristic function C(j3'; s') = C k (j3'; s') of the quan- 
tum state of the initially excited cavity mode and the 
characteristic functions C a {fio\ s a ) = C ka ((3 a ; s a ) of the 
quantum states of the incoming field (<r = in) and the 
dissipative channels (tr = A) as 



c out (/M;s) = e-2^ |2 c[v^M 



(117) where 



Z(t)=r){t)s' + £)Ca(tK-« 



(124) 



(125) 



A = Pi{t) = / dtU&(LU,t)[3{L0). (118) 

Let us further consider the case, when the non- 
monochromatic modes of the incoming field and dissi- 
pative channels corresponding to B^\t), i ^ 1, are in 
vacuum state at the initial time io- Then we may as- 
sume that the resulting characteristic function factorizes 
as well, with 



exp 



ftS£Jt(*)-H.c. 



(119) 



being the characteristic function of the relevant outgoing 
mode. Using Eq. 1(109(1 and noting that the commutation 
relation 



a k (t Q ),B^(t)]=0 (t>t ) 



(120) 



holds (see Appendix |UJ), we may rewrite Eq. I(119|) as 



Cfcout (/?i , t) 



exp 



x exp 



PiB^it) ~ H.c. 



(121) 



Noting that according to Eq. I(lll|) B^^ (t) is a func- 
tional of &fci n (w,<o) and CfcA(w,^o) and assuming that the 
density operator (at the initial time to) factorizes with 
respect to the cavity field, the incoming field, the dissi- 
pative channels, we obtain 



Cfcout t) 



exp 



AB< 1)f (t)-H.c. 



x (exp 



(122) 



Making use of the commutation relations ((70(1 and 



B^(t),B^(t)]=l-r, k (t), 



>(i)t 



(123) 



[^^(iU^NCfc^)]- 
From Eq. ((124() the phase-space function in s order can 

be derived to be 

2 1 

P ont {a,t;s) = --T7-r 

x J d 2 a' P(a';s')Y[J d 2 a a P a ( ary ;s a ) 



x exp 



£(*) 



provided that 



£(*) > o, 



- a 



(126) 



(127) 



where the equality sign must be understood as a lim- 
iting process. To calculate rj(t) [Eq. 1(108(1 ] and ( a (t) 
[Eq. fTlHl ] we make use of Eqs. ®, P ^ -l(Trfl ]l . ifTHT)! . 
and ((113|> . Straightforward calculation yields [T = T k , 

A ± = A k± , T(°) = T< o) , Ai o) = 4°] , = R[ o) , 7ct = lka , 

r = Tfc] 



, ? (t^oo) = ^, 

6a ( t _Oo) = 2Sj»l + | fl (o)|2 



r 

„(°), 



2ni? 



C±(t^co) = 7rad7± 



r 2 



1-4 



(o)|2 



c Ai o) Air(°) 



2n\l 



c A±* A±T(°' ) 

2ml r 



(cav(£ — > OO) 



(o) 

r 2 



(128) 
(129) 

(130) 
(131) 
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where the damping rates 7 ra d, 7abs, and jx are defined 
according to Eqs. J72J) and Eq. J73J), and 



(°) _ (o) 
^rad ^/crad 



: t 

2|mir k 



(o)|2 



(132) 



Eq. ((126(1 together with Eqs. 1(128() - ((131|> generalizes the 
result in Ref. since it fully takes into account the 
noise associated with the dissipative channels. 



C. Thermal noise 

Let us consider the typical case of the dissipative chan- 
nels being in thermal states, i.e., 



W x {a) = - 



1 



-2|a| 2 /(l+2nx) 



7T 1 + 2n A 



(133) 



(fix, average number of thermal quanta) and calculate 
the Wigner function of the quantum state of the relevant 
outgoing mode. Inserting Eq. (|133|) into Eq. ((126(1 . after 
having set sx — therein, performing the ax integrations, 
and setting s = s' = Sj n = 0, we derive 



W out (a,t) 



1 



d a 



x exp 



d 2 /3WV)^ in (/3) 
2\y/^)a' + y/&ffl0- 



where 



^(t) = i-^)-a)+2E fl ^(<). 



(134) 



(135) 



Let us consider Eq. (|134|) together with Eq. I|135|l for 

some typical situations in more detail. 

(i) When the incoming field is in the vacuum state (un- 



used input port, Wi n ((3) — 2tt 1 t 



-2|/3|- 



and the dissipa- 



tive channels — in particular, the coupling mirror — are in 
the vacuum state as well, then almost perfect extraction 
of the quantum state of the cavity mode requires the 
condition 



1 - 77(f) 



> 1 



(136) 



to be satisfied. In other words, on recalling Eq. 1)128(1 . 
the nonradiative cavity-field decay rate must be small 

compared with the radiative one, 7abs/7rad < ^l — a con " 
dition that can be hardly satisfied for a high-Q cavity 
presently 0> GU • How small — depends on the nonclassi- 
cal features of the quantum field to be extracted. Notice, 
this condition can be also obtained by means of quantum 
stochastic approach to the problem |39(. 
(ii) When the input port is unused but the dissipative 
channels are thermally excited, then, as one can easily 



see from Eq. I] 134(1 . the condition to ensure nearly per- 
fect extraction of the quantum state of the cavity field is 



i- v (t)+2j: x fix(x 



> i 



(137) 



The condition Eq. ((137(1 strengthens even more the re- 
quirement of smallness of nonradiative cavity-field decay 
rate compared with the radiative one. Particularly, the 
value of J^x ^aCa(*) should be as small as possible to en- 
sure that the effect of thermal noise effectively does not 
play a role. This is obviously the case when both fix and 
(x(t) are sufficiently small. Needless to say that small 
values of fix require sufficiently low temperatures. For 
cavities with high-quality mirrors with the finesse 

of several hundred thousands, the second, the third and 
the forth terms on the right-hand side of Eq. I(13U|) are 
of an order smaller magnitude than the first term on the 
right-hand side of Eq. I(13U|) . as well as Ccav, Eq. ((131(1 . 
and may be therefore disregarded in the sum ^ A n\C\(t). 
That is to say, in case of unused input port dissipation 
due to absorption in the coupling mirror can be effec- 
tively described by adding appropriate Langevin noise 
forces in Eq.lJSSJ. 

(iii) To describe typical problems on engineering of non- 
classical states of light let us assume, that the dissipative 
channels are again thermally excited and the incoming 
field mode with the mode function xiln( w i £) according to 
Eq. 1(113(1 is prepared in some nonclassical state. Then, 
the quantum state of the outgoing field mode is the one of 
the cavity-mode superposed with the reflected incoming 
field mode as well as the modes of the (thermally excited) 
dissipative channels. The weights of the modes of the in- 
coming field and the cavity-mode field in the resulting 
superposition are defined respectively by the fractions 



Cin(i) 



l-V(t)-Cm(t)+2J2 x fix<:x(t) 

V® 

l-r/-(i)-Cin(t) + 2£ A n A CA(*) 



(138) 
(139) 



Notice, that dropping the absorption in the cou- 
pling mirror, |i^ o) | = l, and Eq {Tg^l reduces to 
(in(t — ► oo) = [1 — r)(t — ► oo)] 2 . The additional noise as- 
sociated with the coupling mirror reduces the fraction of 
the input field in the resulting superposition, and, there- 
fore, represents the absorption of the incoming field mode 
in the coupling mirror. For high- Q ca vities with the fi- 
nesse of several hundred thousands |4a.l49j . the unwanted 
losses in the coupling mirror reduce the weight of the in- 
coming field mode by about 50%. In this way, the quan- 
tum state of the output mode carries additional noise. 

Notice, to obtain the above given results we have as- 
sumed, that the nonmonochromatic modes of the incom- 
ing field and the dissipative channels corresponding to 
(t), i 7^ 1, are initially prepared in the vacuum state. 
In practice, this is not necessarily the case, especially 
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with regard to the dissipative channels associated with 
the coupling mirror, due to the finite number of ther- 
mal quanta and the impossibility to prepare the mode 
of a dissipative channel. As a consequence, additional 
noise is fed into the cavity. Moreover, it is straightfor- 
ward to prove with the use of Eqs. (SI) , ^^j) and l|12()jl . 
that there are necessarily more than one (nonmonochro- 
matic) mode functions of the outgoing field, including 
the one corresponding to t (t), that lie in the rele- 
vant frequency interval, defined by the bandwidth Tk of 
the cavity mode resonance frequency uj^. Therefore, if 
the nonmonochromatic modes of the dissipative channels 
and the incoming field corresponding to B^\t), i ^= 1 
are initially prepared in another than the vacuum state, 
the quantum state of the output field in the relevant fre- 
quency interval is a mixture of modes of the outgoing 
field corresponding to tyS, ui (t) including the relevant one, 

which corresponds to Sj^jL. (t) . The mode analysis for this 
case will be performed in detail in a forthcoming paper. 



VI. SUMMARY AND CONCLUSIONS 

Within the frame of exact quantum electrodynamics 
in causal media we have studied the input-output prob- 
lem of a high-Q cavity. Making use of the representa- 
tion of the quantized electromagnetic field in dispers- 
ing and absorbing planar (dielectric) multilayers as given 
in Ref. |47j . we have considered a one-dimensional cav- 
ity bounded by a perfectly reflecting mirror and a frac- 
tionally transparent mirror, which is responsible for the 
input-output coupling. In order to study the effect of 
unwanted losses such as absorption losses, we have al- 
lowed both the medium inside the cavity and the coupling 
mirror to be absorbing, by attributing to them complex 
permittivities. Moreover, we have assumed that there are 
also active atoms inside the cavity, which are supposed to 
interact with the medium-assisted electromagnetic field 
via electric-dipole coupling. 

We have calculated the electromagnetic field both in- 
side and outside the cavity. It has turned out that in a 
coarse-grained approximation, i.e., on a time scale that is 
large compared with the inverse separation of two neigh- 
boring cavity resonance frequencies, the intracavity field 
may be expressed in terms of standing waves, and bosonic 
operators associated with them can be introduced which 
obey quantum Langevin equations. In this approxima- 
tion, the radiative losses due to the input-output coupling 
and the absorption losses can be regarded as represent- 
ing independent dissipative channels, each giving rise to a 
damping rate and a corresponding Langevin noise force. 
The result shows that the Hamiltonian used in quantum 
noise theories |33| to treat a leaky cavity can be simply 
complemented by bilinear interaction energies between 
the cavity modes and appropriately chosen dissipative 
channels to model unwanted losses such as absorption 
losses. 



However, this intuitive concept fails with respect to the 
operator input-output relations in general. As we have 
shown, the absorption losses attributed to the coupling 
mirror give rise to additional force terms in the input- 
output relations which cannot be simply inferred from 
the above mentioned interaction energies between the 
cavity modes and the dissipative channels introduced to 
model the mirror-assisted absorption. Hence the input- 
output relations obtainable from standard quantum noise 
theories would be incomplete. 

Finally we have used the exact operator input-output 
relations to explicitly calculate the quantum state of the 
outgoing field as a function of time, assuming that the 
quantum state of the cavity field is known at some ini- 
tial time. To be more specific, we have restricted our 
attention to a single cavity mode and assumed that the 
process of quantum state preparation is sufficiently short 
compared with the decay time of the mode under consid- 
eration, so that the time scales of quantum state prepa- 
ration and extraction from the cavity are well separated 
from each other. Introducing the relevant modes of the 
incoming and outgoing fields, i.e., the modes the cavity 
mode couples to, we have expressed the s-parameterized 
phase-space function of the quantum state of the relevant 
outgoing mode in terms of the phase-space functions of 
the quantum states of the cavity mode, the relevant in- 
coming mode, and the dissipative degrees of freedom re- 
sponsible for unwanted losses. It should be mentioned 
that the generalization to more than one cavity mode 
initially excited is straightforward. 
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APPENDIX A: RECURSION FORMULAS FOR 
FRESNEL COEFFICIENTS 

To calculate and tij , we first note, that in the case 
\i — j\ = 1, i.e., single-interface transmission, they are 
defined according to 

r « sr */' = frf = _rj/i ' (A1) 

t i j=ti/j = 0-+r i / j ) = §-t ji , (A2) 
Pi 

leading to 
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In the general case, the relations (see Ref. |5(j) 



Ti/j/k 

t 



1 - ^/i^/fce 2 ^^ 



(A4) 
(A5) 



[min(i, k) < j < max(j, k)] hold. With these formulas at 
hand, we can calculate recursively all the quantities r^/j 



and ti/j, since 



nk = r i/k = r i/j/k , 

tik — ti/k ti/j/k 



(A6) 
(A7) 



for any j with min(«,A:) <j < max(i,fc). Note that in 
the system under consideration r\o=— 1, since perfect 
reflection from the left-side mirror of the cavity field has 
been postulated. 



APPENDIX B: DERIVATION OF EQ. ifffill 

Inserting Eq. \2'2\) in Eq. (|19f) for j = 3, we may write 



E^'(z,u,t) at z = 0+ (cf Fig. 03) as 



p(3) , .s| 
£=Lout ,free l z ) w ) r J | 2= Q+ 



+ 



Di 

t 23 e^ d 
D 2 



where C^±(l, uj, t), C±' (u>, t), and Cf (uj, t) are given by 
Eqs. lO, ED, and respectively. Using Eqs. 
(|33|1 . and (|A4|) . the following relations can be easily 
proved: 



C^ + (l,UJ,t)-C^_(l,UJ,t) 



C { +\uj,t)+r 20 C m (uj,t) + r 30 Cf(uj,t), 

(Bl) 



V(3), 



A(2) 



(3), 



f"30 - ^31 



W 2 t 23 

h 3 t, 
-Dl 



131-31 2i0 x l 



(B2) 
(B3) 



Combining Eq. (|B1|> with Eqs. (|B2f> and (|B3(I . we arrive 
at Eq. 



APPENDIX C: DERIVATION OF EQS. lffl|l - (f?3|l 

We solve Eq. (2} [equivalently, Eqs. J3fJ) and l|3"7])] by 

iteration to obtain T k in leading order as 



2ni 



r(l - |r 13 | 2 ), 



(CI) 



where ni, r\ 3 , and the parameter introduced in the fol- 
lowing are taken at the (unperturbed) frequency uj k . By 
means of 



ri3 



-r 2 i + r 23 e 



2if3 2 d 2 



1 - r 23 r 21 e 2i ^ d ^ 
Eq. IClfl can be rewritten as 
1 4ni 



(C2) 



\D'\ 2 \n 1+ n 2 \ 2 



n' 2 (l-|r 23 |V 4 ^) 



//; . ( r 23 e - 2i /35 d - r 23 e 2 ^ 2d 



(C3) 



Further, from Eqs. (|4H (w = ojfc) and (|A5|) we find 



16ni|n2| 2 n3 



£ 2 | 2 K+n 2 | 2 |n 2 +n 3 



Making use of Eqs. I|47|l and (|48|l . we derive (uj — uj^) 



-p> 2 'd 



D' 2 \n 1 +n 2 



x [„£ (e^ d - e-^ d ) (l + \r 23 \ 2 e~ 2 ^ d ) 
- m' 2 ' (e^ d - e"^ d ) (r 23 e^ d + r^e-W* 
Thus, combining Eqs. (|C3|1 - HC5(1 . we arrive at 



r fc = 



2 ml 



(\T k \ 2 + Y. 



A 



(C5) 



(C6) 



which matches Eq. J7TJ together with Eqs. (|?2"j) and JJJJl. 



APPENDIX D: PROOF OF THE 
COMMUTATION RELATION (17011 

To prove the commutation relation l|7()|) . we recall the 
definition of (z, i), namely 



E { k 1 \z,t)= I duE^'(z,uj,t), (Dl) 
J(A k ) 



where if \z, oj, t) is defined by Eq. I|19|) for j = 1- Using 
Eq. I|21(l , employing the integral relation 



ImG(zi, z 2 , uj) — — j dx e" (x,uj)G(zi, x,uj)G* (z 2 ,x,uj), 

(D2) 

and recalling the commutation relation Q), after some 
algebra we find that 



[E ! i\z 1 ,t),E^\z 2 ,t 



(i)t. 



= Skk' 



f d^^ImG* 11 )^,^^). (D3) 
J(A k ) *A 



To perform the integration we extend 
the lower (upper) integration limit to 
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— oo and rewrite ImG' 11 '^, Z2, u>) 

= [G^(z 1 ,z 2 ,oj)-G^*{z 1> z 2 ,u ) )]/(2i). Then, re- 
calling the definition of G^(zi,z 2 , w) [G (11) *(zi, z 2 , w)] 
from Eqs. 1221) - H26|) with j = 1, we evaluate the integral 
applying the residue theorem for the poles determined 
by the zeroes of the function D\(uS) [D\(u))]. Thus, for 
sufficiently high-Q cavities, Tk <C we obtain 



[E^( Zl ,t),E^M] 



n-k' 



e \ei\lA 



sin^i^fe)^] sin[/3J , (w fc )z 2 ]. (D4) 



Comparing Eq. (|D4fl with Eq. i|66|) [together with 
Eq. I(j7|) ]. we then easily see that the commutation re- 
lation J7UJ) holds. 



APPENDIX E: PROOF OF THE COMMUTATION 
RELATION <f93j) 

From Eq. (H3J| it follows that 

[hout(t),b{, ont (t')] = ±- / do,/" dc/|a out | 2 

2n J(A k ) J(A h ,) 



:U =[^ t (0 + , W ,t),l21(0+,c',i')]. 



(El) 



which in the source-quantity representation reads as [cf. 
Eqs. (HU and (E3>] 



[hout(t),bl out (t')] 
= ±[ to 



to 



/iocSvw 



x { [li 3 out, fr cc(0 + , w, *), liw, frec (0+ a;', 0] 

+ [i£V,M > ii2V,"',* / )] 

+ [lil, free (0+a;,t),iS t (0+a;',t / )] 

+ [li? (0+ i) ) l2 ) o t ut , free (0+ a;', 0] } ■ (E2) 

Using Eq. (|20|l one easily finds 



1 wfwf 



[JL (0 + ,wi,*i),iL (0 + ,W2,*a)] = , , ,, , 
x ^ Im G (13) {z A > , 0+ wi)Im G (13) (z A ,0 + , w 2 ) 



AA' 



x / at' / dt" e(*i - 1')©(*2 - 1") 

x e-^i^-tOe^^-t")^,^^^^')]. ( E3 ) 



Further, from Eqs. l|T5 |l -l(T7 |l it follows that 



[f{ Iee {z',Lll,ti),dA{t')] = -Mo^i 



2 eo 



x J2 I d* /, ©(* / -* // )G*(*A',« , ,wi)e _fa ' l(tl -*" ) 

x [4, (f"),d A (0]- ( E4 ) 

At this stage we first multiply both sides of this equa- 
tion by iu}l(j, y/he /(irA) e"{toi) G (3:,) (0 + , z',u{) and 

perform the sum an d the integrals J*r .-, dz', by 

making use of Eq. (|D2|I . Next we multiply the result by 
-i/(7:e A)(uj/c 2 )e(t 2 - t')lmG^ 3 Hz A ,0+,Lu 2 )e^ (t ^ t '\ 
take the sum with respect to A and the time integral 
with respect to t', and recall the free-field and the 
source-field definitions (|19f) an( l jHH together with 
Eq. leading to 

- ("3") - (3lt 

[^fefrcc (0+ wi,ti),^'(0+ wa.ta)] 

1 ^1^2 



TT^-A 2 c 4 



^ Im G< 31 > (0+ z A , , wi )Im G< 13 > 0+ w 2 ) 



AA' 



x y d^ y dt" e(t a - 1')©(*' - o 



x e 



-iwi(ti-t")pto>2(ta-0 frf * , r+ ""I ^„ (V 



[d A >(t"),d A (t')]. (E5) 



Using Eqs. (|E3f> and l|E5J) we then derive, on recalling 
that e(x) + Q(-x) = l, 

* (3) * (3")t 

+ [ElL(0 + ^i,h),El "(0+ |W2> t a )] 

+ [E® (0+ Wl| ix) , i^Ie ^, t2)] 

1 w 2 w| 



7r 2 £ 2 „4 2 c 4 



^ Im G (13) (z A , 0+ wi)Im G (13) 0+ w 2 ) 



AA' 



x y d*' y at" [e(*i - o©(*' - *")©(*" - 1 2 ) 
+ e(t 2 -t")Q(*"-*0©(*'-*i)] 

x e -^iti-t') e i"^-t")[dA>{t'),dA{t'')\. (E6) 
In a similar way one can calculate the commutator 

- (3) - (31 1 

lE} k > (0+ u, i) , ^L ;T (0+ w', O] ■ where 



=,(3) 



I^ lfnJ e(«,w,t) = e-* 2 Gi 3) (^,t), (E7) 

with C^\uj,t) being given by Eq. (O- That 
is, multiplying both sides of Eq. I|E4|I by 
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{ c/ [2n 3 ( Wl )]} y/he /(irA) *M(ux) e^ 3 ( Wl ) z ' and 
[i/(7reo^)](^/ C 2 )e(t 2 -t')ImG( 13 )(^,0+w 2 )e^^- i ') ) 
taking the integral J, g -, dz', the time integral with respect 
to i', and the sum with respect to A, we arrive at 



[l2iUee(0 + , W l! t 1 ),l^ T (0+ W 2 ,t 2 )] 



(3)t, 



■K 2 elA 2 c 4 2cn 3 



dz'e^wi) 



x Y, e l ^ l)z 'G^ 3 >(z A ,,z, Wl )Im G^ 13) (z A: 0+ c^ 2 ) 

x y dt' /" dt" 0(i a - t')Q{t' - t") 
x e -<"i(*i-t") e *^(ta-*')[d A ,(t") j d il (t')]. (E8) 
Now we may calculate the commutator 

- (3) a (3)t 

[Jw,fr ee ( 0+ > w >*)>£fc s (0 + ,uj' ,t% using the iden- 
tity 

[li^ frcc (0 + ,c,t),^ 3 s )t (0+c',i')]. (E9) 



Combining Eqs. (|E2]> . (fE5|) . I|E6]) . p8jl . and {E2J), we 

derive 



[SfcootW.S^^)] 
1 

~ 2^ 



(A fc ) J(A fc ,) 



do/ |a OIlt | 2 - 



7T„4 



A ^ A (3-) 

'out, free (O+w',0]- (E10) 

It can be shown |47j that the operators Sout.freeO^) t) 

[and &Q Ut f rcc (w, t)] defined according to Eq. with 

A (3) * (3) 

^out froc( z > w > *) m pl ace °f ^o U t( z , obey the Bose 
commutation relation 

[Wcc(^ t), iL^iu/, t')] = e-W-^Siw - J). 

(Ell) 



Hence from Eq. (|E10|) it follows that 



W*),6Lut(*0] =^"^ / dwe-Mt-O. (E12) 



(A k ) 



Extending, within the approximation scheme used, the 
limits of integration to — 00 and 00, we arrive at the com- 
mutation relation (|93J) to be proved. 



APPENDIX F: DERIVATION OF EQ. (|94jl 

To prove Eq. I|94|l , we first show that the integral term 
on the right-hand side of the Eq. I|88|) can be rewritten 



as follows: 



H fc (t)= / dt'e(t-t') 



1 /.t+At 

= J_ / dw' / dt" 

2 "" y A a) t Jta 



,i(w'-u)(t-t") 



1 Au k Jt 

x / dt / e(t-i / )e-* n *(*-*'> 



Tfc(w')6fei n (w', t') + ^ A fc A(w')cfcA(a/, t') 



(Fl) 



(At^> A fc 1 ). To verify, we first perform ^''-integration on 
the right-hand side of this equation to obtain 

E fe (i) = — [ do/ / dt' e(t - t') e - <n * (*-*') 

gi(t»;— u/)At 1 + 1 e i(u-uj')(t<i-t) 

v 

i(w — a/) 

T^h^Lu'^^ + YAkxi^CkxiLu^t') . (F2) 

A 

In the coarse-grained approximation used, i.e., t> Alu^ 1 
(t = At, t — to), we may let 

g<(u)-u/)T _ ^ 



i(w — cj') 



iP- 



£P- 



w — w 
1 



7 + ttS(lo - u') if r > 0, 



(F3) 



w — UJ 



- - ttS(lo — u') if r < 



(P, principal value). Now the w'-integration can be eas- 
ily performed to see that Eq. (|F1|) is correct within the 
approximation scheme used 

Making [on the right-hand side of the Eq. I|F1|| ] the 
change of variables: t' — > t — t" + t' and performing u>'- 
integration [T k — T k (uj k ), A kx = A kx (uj k )}, we find 



Sfc(t) 



'2tt 



t+At 



dt" / dt'e(t"-t') 



x e -iU k {t"—t') e —iu(t—t") 



x T k b kin (t')+J2 A kxc k x{t' 

X 

Comparing Eq.l|F4|l with Ea.l(BH)l. we see that 
S fc (t) = 



(F4) 



2ni(u k )l 



'2tt 



t+At 



dfe-^-^a^t"). (F5) 
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Substitution of Eq. (|F5p into Eq. eventually yields 
Eq. JM|). 



APPENDIX G: DERIVATION OF EQ. Jl20l) 

Using Eq. ((SSJl together with Eqs. and the 

commutation relations H50JI and l|51|) it is not difficult to 
prove that 



[a k (t),b{, in (u,t')] = Skk' 
[ak{t),cl, x (uj,t')] = S kk , 



_2n 1 (u k )l 



2ni(uj k )l 



* T k i e -Mt-f) 

^/2~n u — fife 

(Gl) 



^/2tt u — Ofe 



Using Eq. I|110() together with Eqs. and (|107|l . from 
Eqs. i|Gl|) and (|G2[I we can calculate the commutator 
(I12()|l in a straightforward manner. 
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